Higher-dimensional puncture initial data 
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We calculate puncture initial data, corresponding to single and binary black holes with linear 
momenta, which solve the constraint equations of D dimensional vacuum gravity. The data are 
generated by a modification of the pseudo-spectral code presented in [l[| and made available as 
the TwoPuNCTURES thorn inside the Cactus computational toolkit. As examples, we exhibit 
convergence plots, the violation of the Hamiltonian constraint as well as the initial data for D — 
4, 5, 6, 7. These initial data are the starting point to perform high-energy collisions of black holes in 
D dimensions. 



I. INTRODUCTION 

Numerical relativity in higher-dimensional spacetimes 
could be a powerful tool to study a variety of physical 
concepts, such as the stability of black hole solutions and 
their interactions, as well as for producing phenomeno- 
logical information of relevance for TeV gravity scenar- 
ios f2U3|. In such models, the fundamental Planck scale 
could be as low as 1 TeV. Thus, high-energy colliders, 
such as the Large Hadron Collider (LHC), may directly 
probe strongly coupled gravitational physics |5l-[l0l|. In- 
deed, in this scenario, particle collisions could produce 
black holes 0, Moreover, the production of black 
holes at trans-Planckian collision energies (compared to 
the fundamental Planck scale) should be well-described 
by using classical general relativity extended to D di- 
mensions (see [111] and references therein). Numerical 
simulations of high-energy black hole collisions in higher- 
dimensional spacetimes, then, could give an accurate es- 
timate of the fractions of the collision energy and an- 
gular momentum that are lost in the higher-dimensional 
space by emission of gravitational waves; such informa- 
tion would be extremely important to improve the mod- 
eling of microscopic black hole production, and of the en- 
suing evaporation phase, which might be observed during 
LHC collisions. 

The growing interest in dynamical aspects of higher- 
dimensional spacetimes led to the development of higher- 
dimensional numerical relativity [l2l - [l5j and to the pro- 
duction of the first black hole collisions in higher- 
dimensional spacetimes starting from rest [l6l [T7| , and. 



most recently, of black holes with initial boost p3|j- 
Aside from their immediate relevance in the context of 
TeV gravity scenarios, high-energy collisions of black 
holes provide fertile ground for probing strong-field ef- 
fects of general relativity, such as cosmic censorship, lu- 
minosity limits, hoop type conjectures and zoom-whirl 
behaviour p|-|^. 

The modeling of generic black hole or neutron-star 
spacetimes in the framework of Einstein's theory of gen- 
eral relativity requires a numerical treatment because of 
the enormous complexity of the equations in the absence 
of high degrees of symmetry. For such numerical mod- 
eling, the Einstein equations are cast as a time evolu- 
tion or initial value problem as originally formulated by 
Arnowitt-Deser-Misner (ADM) [2^ and reformulated by 
York [23|. Numerically generating a solution then con- 
sists of two basic steps: (i) the construction of initial data 
which satisfy the constraint equations and represent a re- 
alistic snapshot of the physical system under considera- 
tion and (ii) the evolution in time of these initial data. 
In this work we will focus on the first step and provide a 
formalism for constructing initial data for black hole bi- 
naries with nonzero boosts in higher, D > 5, dimensional 
spacetimes. 

Most work on the generation of initial data in 3 -I- 1- 
dimensional general relativity is based on the York- 
Lichnerowicz split [25l - [28j which rearranges the degrees 
of freedom contained in the three-metric and extrinsic 
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^ These simulations start from superposed single boosted black 
hole data without application of a constraint-solving procedure; 
our work is motivated in part by providing initial data for a com- 
parison and thus calibrating the impact of constraint violations. 
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curvature Kij via a conformal transformation and a split 
of the curvature into trace and traceless part 

1 (1-1) 

followed by a transverse-traceless decomposition of ei- 
ther the traceless extrinsic curvature Aij or a confor- 
mally rescaled version Aij thereof. For details of these 
physical or conformal transverse-traceless decompositions 
as well as the alternative conformal thin- sandwich for- 
malism, we refer the reader to Cook's review |22], Alcu- 
bierre's book [30] and references therein. 

One of the main advantages achieved with this de- 
composition is the decoupling of the momentum from 
the Hamiltonian constraint under the additional simpli- 
fying assumptions of conformal flatness, jij = 6ij and 
a constant trace of the extrinsic curvature K = const. 
Quite remarkably, the resulting equations for the mo- 
mentum constraints admit analytic solutions describing 
multiple black holes with nonvanishing spins and linear 
momenta (sll . There then remains a single elliptic differ- 
ential equation, the Hamiltonian constraint, for the con- 
formal factor which requires a numerical treatment. 
By applying a compactification to the internal asymp- 
totically flat region, Brandt & Briigmann ^32*] derived a 
method of solving the Hamiltonian constraint which is 
now generally referred to as the puncture method. These 
puncture data form the starting point for the major- 
ity of numerical simulations using the moving-puncture 
method [H, [sJl . We note that an alternative approach 
to evolving the Einstein equations based on the gener- 
alized harmonic formulation has been implemented with 
similar success [H, . 

As shown by Yoshino et al. [s^l , the existence of an an- 
alytic solution of the Bowen-York type for the momentum 
constraints carries over to the theory of general relativ- 
ity in higher dimensions. These authors also solved the 
Hamiltonian constraint using a finite difference method. 
Here, we will present a generalization of the spectral 
solver by Ansorg et al. used inside the CACTUS 
toolkit [3^, , that solves the Hamiltonian constraint for 
black hole binaries in D > 5 dimensions with nonvanish- 
ing initial boost, and preserves the spectral convergence 
properties observed in four dimensions. The compactifi- 
cation implemented in this solver further facilitates direct 
interpolation of the initial data onto computational grids 
of arbitrary size with mesh refinement employed in state- 
of-the-art numerical simulations of black hole spacetimes. 
Finally, the initial data are generated in variables that 
can be straightforwardly translated to evolution systems 
common in numerical relativity, such as the Baumgarte- 
Shapiro-Shibata-Nakamura [io, |4l| system in Cartesian 
form. In our case, we specifically translate the spec- 
tral solution into the formalism of which can thus 
be readily used to perform high-energy collisions of two 
black holes in D dimensions. We complement our study 



with the analytic perturbative analysis of boosted single 
puncture initial data. 

This paper is organized as follows. In Sec. [Hi we re- 
view the constraint equations, the Brill-Lindquist [l^ 
and Bowen-York [sij type initial data and introduce an 
appropriate coordinate system. In Sec. Illli we provide 
the explicit form of the elliptic equation that must be 
solved for determining the Bowen-York initial data that 
describes a boosted head-on collision in D dimensions. 
We also present its dimensional reduction to four space- 
time dimensions, following In Sec. llVi the case of 
a single black hole with linear momentum is considered 
and it is shown that an approximate analytic solution 
can be found for all variables, including the conformal 
factor, in the limit of small momentum to mass ratio. 
In Sec. |Vl we explain the modifications we have coded 
to study the case of two black holes with aligned linear 
momentum P/r^'^ = ±0.8 (rg is the Schwarzschild ra- 
dius) and present results for convergence tests, constraint 
violation and the solution itself in D = 4, 5, 6, 7. 

Notation: In the remainder of this work, early lower 
case Latin indices a, b, c, . . . extend from 1 to D — 1, 
late lower case Latin indices i, j, k, . . . run from 1 to 3 
and early upper case Latin indices A, _B, C, . . . from 4 to 
D-l. 



II. HIGHER-DIMENSIONAL INITIAL DATA 

The starting point for our discussion is a {D — 1)- 
dimensional spacelike hypersurface S with induced met- 
ric 7ab, and extrinsic curvature Kah embedded in a D- 
dimensional spacetime. By generalizing the ADM de- 
composition, the spacetime metric is given in the form 

ds^ = -a^dt^ +^ab{dx'' + P''dt){dx^ + (3^dt) . (2.1) 

We conformally decompose the spatial metric and extrin- 
sic curvature as 

_ 4 

lab = ■0°"^7ab , 

1 _ - (2-2) 

Kab = V' '^Aab + jy—jJabK , 

which generalizes the 3-1-1 dimensional Eq. (|l.ip . Note 
that we here represent the traceless part of the extrinsic 
curvature by the conformally rescaled version A"''. We 
assume that the conformal metric jab is hat and impose 
the maximal slicing condition, K = 0, which lead to the 
decoupling of the constraints mentioned above. With 
this choice, the higher-dimensional initial data equations 
in vacuum become [s^, Ell 

aai"" = , (2.3) 

A^+ ^^^~^^^ ^-^i°''i.6 = 0, (2.4) 

where i"'' ee ^""■^'"^Acd, A = dad" is the flat space 
Laplace operator and the first equation holds in a (£> — 1)- 
dimensional Cartesian frame X"^ = (x^, x^, . . . , x^~^). 
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A. Brill-Lindquist initial data 



M(i) of the z*'' hole by 



For the time-symmetric case Kat = 0, Eq. (I2.3P is au- 
tomatically satisfied, and Eq. (j2.4p reduces to the D — 1- 
dimensional flat space Laplace equation, 



Ai/' = 



(2.5) 



For asymptotically flat spacetimes, the confornial factor 
satisfies the boundary condition 



lim -0 = 1 



and a solution to Eq. (|2.5|) is given by 



(2.6) 



N 



Hi) 



i=l (i) 



(2.7) 



where r^j) = \r — a;(i)|, x^j) is the (arbitrary) coordinate 
location of the z*^ puncture, and the mass parameter 
is related to the horizon radius rs^^^ and the ADM mass 



Hi) 



„-D-3 



16nM, 



Ad-2{D~2) 



here Ad-2 is the area of the unit [D — 2)-sphere and we 
have set the D dimensional Newton constant to unity. 

These closed-form analytic data are the D-dimcnsional 
generalization of Brill-Lindquist data |42| and describe a 
spacetime containing multiple nonspinning black holes at 
the moment of time symmetry, i. e. with vanishing linear 
momentum. 



B. Bowen-York initial data 

In order to numerically evolve black holes with nonzero 
boost, we need to generalize Brill-Lindquist data to the 
non time-symmetric case. In four dimensions, this gen- 
eralization is given by the Bowen-York extrinsic curva- 
ture, a nontrivial analytic solution of the momentum con- 
straint Eq. (|2.3p . As shown by Yoshino et al. [s^, we 
can write a solution of Eq. (|2.3p describing a spacetime 
of arbitrary dimensionality D containing N black holes 
in the form 



N 

1=1 



(2., 



where 



Aab _ 



At:{D-1) 1 



[D - 2)A 



D-2 r 



D- 



(2.9) 



Here we have introduced = — ^"^^ param- 
eter Pj°j corresponds to the ADM momentum of the z*^ 
black hole in the limit of large separation from all other 
holes. 

In order to obtain a complete set of initial data we 
still need to solve the Hamiltonian constraint (|2.3|) with 
Aab given by (j2.8l) . For this purpose, we follow the stan- 
dard decomposition of the conformal factor into a Brill- 
Lindquist contribution i/jbl given by (|2.7p plus a regular 
correction u 

V' = V'BL+u. (2.10) 
Equation (|2.4p then takes the form 

hu+ ~ A^^Aabir^ . (2.11) 

4(1) -2) ^ ^ ^ 

As in D — A, the higher-dimensional extension of Bowen- 
York extrinsic curvature data can also accommodate an- 
gular momentum of the black holes. In the present work. 



however, we shall focus on initial data for nonspinning, 
boosted black holes only. 



C. Coordinate transformation 

In summary, the initial data are determined by (i) the 
extrinsic curvature Kab obtained by inserting Eq. (j2.9p 
into (|2.8p and the resulting Aab into Eq. (12. 2p . and (ii) 
the spatial D — \ metric ^ab obtained by numerically 
solving Eq. (|2.1ip for u which gives the conformal factor 
via Eq. ((?1U1) and the metric through Eq. ((^ . 

For the numerical solution of Eq. (|2.1ip . it is conve- 
nient to transform to a coordinate system adapted to the 
generalized axial symmetry SO{D — 2) in D = 5 dimen- 
sions and SO{D — 3) in D > 6 dimensions as discussed 
in Sec. I C of Ref. [IJ. For this purpose we consider the 
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(flat) conformal metric in cylindrical coordinates 

^abdx°-dx^ = dz^ + dp^ + {d^p^ + sin^ LpdVLo-i) , 

(2.12) 

where dfln-i is the metric on the {D — 4)-sphere. Ob- 
serve that (/? is a polar rather than an azimuthal coordi- 
nate, i.e. (fi € [0,n]. Next, following we introduce 
"incomplete" Cartesian coordinates as 



p cos ip . 



y = psvaip , 



(2.13) 



where —oo < x < -\-oo and < y < -l-oo. The D- 
dimensional initial data for the spatial metric is then 

^abdx^dx^ = ipT^ [dx'^ + dy'^ + dz^ + y^dVLD-i\ ■ 

(2.14) 

We can transform the D — 1-dimensional Cartesian co- 
ordinates X"^ = {x^, . . . , x^~^) to the coordinate system 
y"" — {x, y, z, ^1, ^2, ■ • ■ , C-D-4) with hyperspherical coor- 
dinates Ci, • • ■ ,C-D-4 by 

x^ — X 

x^ — y cosfi 

x^ = z 

x'* = y sin^i cos^2 {D > 6) 

a;^ = ysin^i sin^2 cos^3 {D>7) 



„D-3 



„D-2 



„D-1 



y sin ^1 • • • sin ,^_D_.6 cos ,^_d_5 
y sin ^1 • • • sin ,^i5_5 cos S^o-i 
ysin^i • • • sin,^_D_4 



iD>7) 

{D > 6) 

{D>5) 

(2.15) 

Without loss of generality, we can always choose co- 
ordinates such that the black holes are initially located 
on the z axis at zi and Z2 and have momenta of equal 

I 



magnitude in opposite directions P^°^ = —P{2)- Inserting 
the momenta into Eq. (|2.9I) then provides the conformal 
traceless extrinsic cuvature and the differential Eq. (j2.1ip 
which is solved numerically for u. 

The class of symmetries covered by the formalism de- 
veloped in Ref. [l3| includes head-on and grazing colli- 
sions of nonspinning black holes with initial position and 
momenta 



x^i) = (0, 0, zi, 0, . . . , 0) , x^2) = (0> 0, Z2, 0, . . 

= (p^o,p^o,...,o) = 



r>a 



.,0) 
(2.16) 



Note that a nonzero is not compatible with the as- 
sumed symmetries. On the other hand, the x-axis can al- 
ways be oriented such that the collision takes place in the 
xz-plane. Our formalism therefore covers general grazing 
collisions of nonspinning black hole binaries in D dimen- 
sions. 

III. FOUR DIMENSIONAL INITIAL DATA FOR 
A GENERAL D HEAD-ON COLLISION 

For illustration and numerical testing, we will in the 
rest of this paper discuss in full detail the case of black 
holes with momenta in the z direction, that is, the case 
given by setting P^ = in Eq. (|2.16p . The linear mo- 
menta are thus given by 



P(i) = (o,o,p^o,...,o) 



pa 



(3.1) 



The rescaled trace-free part of the extrinsic curvature for 
such a configuration is 



1(1) 



A 



(2) 



ab ' 



(3.2) 



where A^^^ and A'f^ are given by Eq. (1^ with ([^1^ and 
p.ip . Using Eq. (j2.15p we can write this in the coordinate 
system y"" adapted to the spacetime symmetry: 



with 



4(1) 

^ab 




Att{D - 1)P^ 









-2)AD-2{x^+y^ + {z-z^Y)^ 




"■AB 






zi)^](z~zi) {D-3)xy{z-zi) 


X 





(D-3)xy{z~zi) 
x[x^+y^ + (D-2){z-zif 



-[x''-(D-4:)y^ + (z~zif]{z-zi) y[x^+y'' + {D-2)(z-zif] 
y[x^+y^ + {D-2){z-zif] [x^ +y^ + {D-2){z- ziy^]{z- zi 



(3.3) 



(3.4) 



and 

4s = - ^1) [2^' + y^ + {z- zi)2] hAB , (3.5) 

where Hab is the metric on the {D — 4)-sphere. The 
expression for A^j, is analogous, but with Z2 in place of 
zi and — P^ in place of P^ in Eq. (|3.3p . 



The formalism developed in ^14] for P)-dimensional 
spacetimes with SO{D — 2) or SO{D — 3) isometrics de- 
scribes the spacetime in terms of the traditional three- 
dimensional metric jij and extrinsic curvature Kij cou- 
pled to a scalar field A and its conjugate momentum K\; 
cf. Eqs. (2.14), (2.26) in [ll|. These are the variables 
evolved in time and therefore the variables we ultimately 
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wish to construct from the initial data calculation. For 
their extraction we first note that 7^ , Kij and Kx are re- 
lated to the {D— l)-diniensional metric ^ab and extrinsic 
curvature Kab by 



lij = lij , lAB — ^hAB , 

7iA = , 



(3.6) 



D - 4:K\ 

K,A=0, K = K+ — ^ . (3.7) 



Using these relations and Eq. (2.14) of [ij] we can express 
all "3+1" variables in terms of those describing the initial 
data 



{D - A)Kx 



K,, = ir\^ + i(f ) , Kx = 2r^y^{P+ + P-) , 



K 



where 



2A 



4tt{D- l)P^(z-zi) 



(D - 2)AD-2ix^ + y^ + {z~ zi)2)^ ' 



iTT{D- l)P^(z-^2) 



P+ 

The conformal factor is 

Ml 



(3.8) 



(3.9) 



^ = 1 + 



4[x2+?;2 + (2_2^)2]P-3)/2 



Ai2 



(3.10) 



2](D-3)/2 



^x"^ + y"^ + {z - Z2) 



and u is the solution of the equation 



(3.11) 

I 



where 



-{D - A){P+ + p-f 



(3.12) 



Our numerical construction of the function u will be 
based on the spectral solver developed in [l[ . This solver 
employs coordinates specifically adapted to the asymp- 
totic behaviour of u at spatial infinity. In order to in- 
vestigate this behaviour, we next consider a single black 
hole with non-zero linear momentum. 

IV. SINGLE PUNCTURE WITH LINEAR 
MOMENTUM 

For a single puncture with momentum P^ located at 
the origin z = 0, Eq. (|2.9p implies 



167r2(D- 1) 



-Pt 



(£» - 2)2yl2^_2r2(-D-2 

so that Eq. (|3.1ip takes the form 



2 + D{D-3) 



87r2(i?-l)2(Z)-3) „2 



{D - 2)M25_2r2(-D-2) 



Pfx 



(4.1) 



(4.2) 



r 



It turns out to be convenient for solving this differential 
equation to introduce a hyperspherical coordinate system 
on the D — 1-dimensional spatial slices, such that the flat 
conformal metric is 

^^2 = ^abdx'^dx^ 



with cosi5 = 
nate 



We further introduce the radial coordi- 



X 



47-D-3 



(4.3) 



which reduces to the coordinate A of Eq. (31) in [l[ for 
the case of _D = 4 spacetime dimensions. Expressed in 
the new coordinate system, Eq. (|4.2p becomes 



dxx + x^^'^ {D- 3)2^2(1 - xy 



P. 



-a(-^j X^— (l + uX) T^(^l+ 2 



d^^ + (Z? — 3) cot ddi) 
D{D - 3) 



1 



sin^ {} 



{d^^ + {D - A) cot Lpd^) 



(4.4) 



with 



1287r2(i:)- 1) 



For D = 4, we recover Eq. (40) of [1]. In order to study 
the behavior of the solution at spatial infinity, we now 



iD-3){D-2rAl_ 
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perform a Taylor expansion in v = ^ 



2i, 



(4.5) 



Odd powers of v have to vanish in order to satisfy 
Eq. (|4.4I) . We have the fohowing equation for ui 



dxx + ^dx + _ 3)2X2(1 - xy 



[di)i) + (£> - 3) cot ddi)] \ u 



D{D-3) 



(4.6) 



In order to solve Eq. (j4.6p . we make the ansatz 



ui = f{X)+g[X)QD{cosd) , 



(4.7) 



where Qoicosd) = {D — 1) cos^ — 1. By solving 
Eq. (j4.6l) . we find that the functions f{X) and g{X) take 
the form 



f(X) 



327r^(Z)~3) 



1 _ XT^ 



(4.8) 



r 



5(^) = 



ki 



X 



1-X 



+ ^2 I — 7^1 - a^rTTTT-TTTTT— rrX 



'2(£)+l)(D-l) ■ 



p+1 

X— 



D— 1 2 2 

^ '(1-X)— 



_D-3' Z5-3' -^15-3'^ 



J^X^il-X)^^,FJT^,^;3^;X 



(4.9) 

where 2^i(a, &;c;X) is the hypergeometric function and 
ki^2 are constants to be fixed by imposing that g(X = 
1) = and g{X = 0) is smooth. Requiring analyticity 
at X = and using the property F{a,b,c,0) — 1, we 
immediately find k2 = 0. 

We are now interested in the large X — > 1 limit. There- 
fore, we use the z ^ 1 — z transformation law for the 
hypergeometric functions pi ]. 



F{a-c+l,b~c+l,2~c, z) ^ (l-z) 



T{2-c)T{a + b-c) 



r(a- c + l)r(5- c + 1) 
r(2-c)r(c-a-6) 

r(i-a)r(i-6) 



F{l-a,l-b,c~a~b+l, 1-z) 
F{a~c+l,b-c+l,-c+a+b+l,l-z) . (4.10) 



Requiring a regular solution, we find that fci has to satisfy 



• D ^4 



ki 



Mtt^D{D ~ 3) 



2 r 



2{D-2] 
D-3 



(D-2nD + l)Al_2 



(4.11) 



Let us write these functions explicitly for D = 4, 5, 7 
(for D = 6 the hypergeometric function does not sim- 
phfy): 



fix) 

9{X) 



^l^[84(l-X)log(l-X) 
+ 84X - 42X2 _ -^4^3 _ rjj^4 

- 4X^ - 2X'^ 



(4.12) 



(4.13) 



These are Eqs. (42-44) in fl*!, with appropriate re- 
definitions. 
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• D 



5 : 

f{X) 
9iX) 



1 - 



16 
9^ 
80(1 



8l7r2X2 



41og(l-X) 



/(^) 



128(l-X2) 



257r4 



(4.14) 



(4.15) 



(4.16) 



28 



, -30J{l-X)X 

1257rV(l - ^)^''' 

+ 40v/(l-^)^^ - 16v/(l-^)^^ 

+ SttX^ + 6 (5 - lOX + 4^2) arcsinVx 

(4.17) 

Analyzing these expressions, we can anticipate the con- 
vergence properties of the numerical solutions obtained 
in terms of pseudo-spectral methods. For instance, ana- 
lyticity of / and g suggests exponential convergence. As 
will become clear in the next section, we are interested 
in the convergence properties in a coordinate A behaving 
as A ~ 1 — i, for large r. We thus introduce a coordinate 
A that satisfies 



X- (l + (A-i- 1)^-3) 



(4.18) 



In terms of the A coordinate, we find that the functions 
/ are analytical. For the function g in the vicinity of 
A ~ 1, the leading terms behave as follows: 



• D = 5 

giA) ^ - 

• £1 = 6 

• D^7 



80 
8l7r^ 



(l~A)4[81og(l-A) + 7] , (4.19) 



,(A)^^^(l-A)' 
^^^^ 6272^2(1 ^) ' 



257T- 



.(1~AY 



(4.20) 



(4.21) 



From the behaviour of the functions / and g and Eq. (|4.7p 
we conclude that the first term in the expansion (j4.5p has 
a leading-order behaviour ui ~ l/r^~^ as r — > oo. Iter- 
atively solving Eq. (14. 4p for higher powers of v is compli- 
cated by the presence of the source terms on the right- 
hand side, but under simplifying assumptions indicates 
that higher-order terms Uj > 2 acquire additional factors 
of 1 /r and therefore the leading-order falloff behaviour is 
given correctly by that of ui. This result is confirmed by 
our numerical investigation using finite boost parameters 
as we shall discuss in the next section. 



With regard to the analyticity of the solutions and the 
resulting expectations for the convergence properties of a 
spectral algorithm, we summarize the results of our ana- 
lytical study of a single puncture as follows. In Z5 = 6, 7, 
the leading terms are analytic functions in the vicinity 
of A = 1. Actually, for D = 7, g{A) is analytic in the 
vicinity of any point. Therefore, we expect exponential 
convergence of the pseudo-spectral code. For D = 5, one 
observes the presence of a logarithmic term. This type 
of term is known to arise in JD = 4, when punctures have 
nonvanishing momenta [45l . [4^ and in that case their 
presence makes the convergence algebraic in the single 
puncture case. In the next section we shall investigate 
the impact of the logarithmic terms on the convergence 
properties of our spectral solver. 



TWO PUNCTURES WITH LINEAR 
MOMENTUM 



A. Code changes 

We first explicitly list the modifications applied to 
the spectral solver of Ref. [l| and demonstrate how 
these modifications enable us to generate initial data for 
boosted black hole binaries with convergence properties 
and levels of constraint violation similar to the D = A 
case. For this purpose we start by recalling that the 
spectral solver of [i| employs coordinates 

^e[0,l], Be [-1,1], 0e[O,27r], (5.1) 

which are defined by Eq. (62) of Q, 



X — b 
y = b 
z = b 



2A 1 


-S2 


1-A^l 


+ 52 


2A 1 


-B2 


1-^21 


-hB2 


A^ + 1 


2B 



^2 _ 1 1 + 52 ' 



COS( 



(5.2) 



where b is half of the coordinate distance between the 
punctures. In particular, the coordinate A satisfies 



<^ A^l 



(5.3) 



The first modification consists in adapting the source 
term and Laplace operator according to p. lip . 

Next, we note that the type of high-energy collisions 
which form the main motivation for this work often 
start from relatively large initial separations of the holes, 
\zi — Z2\ ^ rs- In order to obtain high-precision solu- 
tions for such binary configurations, we found it crucial 
to introduce a coordinate A' defined as 



A 



sinh [kIA' + l)/2] 
sinh K, 



(5.4) 



where k is an adjustable free parameter. Note that for 
K = we obtain A = ^{A' + 1) For k > 0, however. 



8 



D 


b/rs 




-3 D — 3 1 D — 3 
Ta Ta 




4 


30.185 


0.8 


3.555 


1.78 


5 


30.185 


0.8 


1.931 


2.27 


6 


30.185 


0.8 


1.415 


2.96 


7 


30.185 


0.8 


1.236 


3.81 



TABLE I. ADM mass obtained with Eq. JS^ in units of 
the "bare" Schwarzschild radius r?"^ = r?"^ + r?~^. The 
variation of the ADM mass with resolution is of the order of 
'^^ for all D and n > 100 grid points indicating that the 



10" 



accuracy in the ADM mass is limited by round-off errors. 



the new coordinate A' provides the spectral method with 
enhanced resolution near ^ ^ 0. 

A further modification is related to the asymptotic 
falloff of the function u as obtained in the previous sec- 
tion, 



„D-3 



(5.5) 



To naturally accommodate this behaviour with the spec- 
tral coordinates used in the code, we have changed the 
variable U of Eq. (5) in [1] to 



{A'-l) 



D-3 



U 



(5.6) 



Note that this U variable is the variable that the code 
actually solves for. 

Finally, we adjust the calculation of the ADM mass 
from the numerical solution. For this purpose, we note 
that, asymptotically 



4r; 



D-3 



4r 



D-3 



1 



4r 



D-3 



„D-3 



Sfer and /.± 



with a = I q — -i rrr^ 

ADM mass is then obtained from 



r^-^ = r'^-^ + r^-^ + 4 lim r^-^u 



(5.7) 
The 



^D-3 



„D-3 



4 -2b 



tanh/ 



D-3 



U{A' - 1) , 
(5.8) 



where we have used Eq. (62) of [1], and Eq. (15. 4p 
and (j5.6l) . We show in Table U the values obtained for 
the ADM mass of some cases we considered. 



B. Results 



We first address the convergence properties of the nu- 
merical algorithm by evaluating the quantity 



Sn,miu) = max |1 - U„/Urn\ 



(5.9) 



where the maximum is obtained along the collision axis, 
i.e. z-axis in our case. Here, the index m refers to a 
reference solution obtained using a large number m of 
grid points while n denotes test solutions using a coarser 
resolution, n < m. The result obtained for black hole 
binaries with initial separation b/rs — 30.185 and boost 
P^/rg~^ = 0.8 in Z? = 4, 5, 6 and 7 dimensions is dis- 
played in Fig. [1] We note from this figure, that achieving 
a given target accuracy Sn.m requires a larger number 
of points n as D increases. We emphasize in this con- 
text, however, that this increase in computational cost 
in higher dimensions is unlikely to significantly affect the 
total computational cost of the simulations which typi- 
cally are dominated by the time evolution rather than the 
initial data calculation. Most importantly, we observe 
exponential convergence up to a level of (5„,m(u) « 10~^ 
for all values of the spacetime dimensionality D. Below 
that level, the two leftmost curves in Fig. [TJ correspond- 
ing to I? = 4 and D = 5, respectively, show that the 
rate of convergence decreases indicating that the logarith- 
mic terms become significant and reduce the convergence 
to algebraic level similar to the observation in Fig. 4 of 
Ref. For D — 6, the convergence remains exponen- 
tial, in agreement with the absence of logarithmic terms 
in the analysis of Sec. IIVI Irrespective of a change to 
algebraic convergence, however, our algorithm is capable 
of reducing the quantity 6m,niu) for all values of D to a 
level comparable to the case D ~ A and, thus, produc- 
ing initial data of similar quality as in 3-1-1 dimensions, 
provided we use a sufficiently high resolution n. 
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We now study the numerical results as obtained for 
D — 4, 5_J), 7 with these adaptations of the spectral 
solver of Throughout the remainder of this section 
we will graphically present results in units of the "bare" 
Schwarzschild radius defined asr^^ = ra" -l-ro". 



FIG. 1. Convergence plot for the b/rs = 30.185, P/r°~^ = 
±0.80 cases. 

For illustration, we plot in Fig. [2] the function u ob- 
tained for the case of 6/rs = 30.185, PVrf-^ = 0.8. The 
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behaviour is qualitatively similar for all values of D, but 
the figure demonstrates the faster falloff for larger D as 
predicted by (|5.5|) . For this plot we have used ua = 300, 
riB — 300 and = 4 grid points. The inset in the fig- 
ure shows the function u in the immediate vicinity of the 
puncture. While the profile develops multiple extrema 
for D > 4, the profile remains smooth for all values of D. 

Finally, we show in Fig. |3l the Hamiltonian constraint 
corresponding to the solutions presented in Fig. [5] as mea- 
sured by a fourth-order finite-differencing scheme of the 
evolution code We emphasize that the violation of 
Eq. (j2.1ip inside the spectral initial data solver is < 10^^^ 
by construction. The independent evaluation of the con- 
straint violation in the evolution code serves two pur- 
poses. First, it checks that the differential Eq. (j2.1ip 
solved by the spectral method corresponds to the Hamil- 
tonian constraint formulated in ADM variables; an error 
in coding up the differential Eq. (|2.1ip could still result in 
a solution for u of the spectral solver, but would manifest 
itself in significantly larger violations in Fig. [3] Second, it 
demonstrates that the remaining numerical error is dom- 
inated by the time evolution instead of the initial solver. 
Note in this context that the relatively large violations 
of order unity near the puncture location in Fig. [3] are an 
artifact of the fourth-order discretization in the diagnos- 
tics of the evolution code and are typical for evolutions 
of the moving-puncture type; see e. g. the right panel in 
Fig. (8) in Brown et al. |48| . 

The solid (blue) curve obtained for the "standard" 
D — A case serves as reference. For all values of D, 
the constraint violations are maximal at the puncture lo- 
cation zi/rg ss 15 and rapidly decrease away from the 
puncture. As expected from the higher falloff rate of 
the grid functions for larger D, the constraints also drop 
faster for higher-dimensionality of the spacetime. 

P/rf--^ =0.8, b/r-s =30.185, n = 300 
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FIG. 2. u function for D = 4, . . . , 7 plotted along the z-axis, 
in units of rg. We used ua ~ ub = n = 300, = 4. We also 
show a zoom around the puncture. 



P/r^-'-' =0.8, b/rs =30.185, n=300 




FIG. 3. Violation of the Hamiltonian constraint along the 
z-axis, evaluated with a fourth-order finite-difference scheme. 
The growth of the constraint violation near the puncture is 
an artifact of finite-differencing across the puncture; see text 
for details. 



VI. CONCLUSIONS 

In this paper we have presented numerical solutions of 
the Einstein constraint equations for the construction of 
initial data containing single or binary black holes with 
nonvanishing linear momentum in D > 4 dimensional 
spacetimes. For this purpose we have modified the spec- 
tral solver of Ref. [jj]. As in Z) = 4 dimensions, the mo- 
mentum constraints decouple from the Hamiltonian con- 
straint under the assumption of conformal flatness and 
a spatially constant trace of the extrinsic curvature and 
allow for an analytic solution describing multiple black 
holes with nonzero momenta. One thus arrives at a sin- 
gle elliptic differential equation for the conformal factor 
or, to be more specific, a regular correction function u to 
the Brill-Lindquist part of the conformal factor. We have 
studied the resulting differential equation in the limit of 
a single black hole with small boost in order to inves- 
tigate the asymptotic behaviour of u at spatial infinity, 
where we find u - For D = 6,7 we further ob- 

serve that u is an analytic function expanded in terms of 
the coordinate A [cf. Eq. (|4.18l) ] around spatial infinity, 
so that a spectral algorithm should provide exponential 
convergence. For D = 5 the expansion of u includes a log- 
arithmic term, but, as has also been observed in Ref. Ij 
for L* = 4, this term is subdominant in the case of a black 
hole binary with equal and opposite momenta. We have 
used the asymptotic behaviour of u to adapt the set of 
coordinates employed in the spectral solver to arbitrary 
dimensionality and further performed a transformation 
to a radial coordinate that ensures sufficient resolution 
near spatial infinity for the case of large separation of 
the two holes. 

The resulting code has been used to calculate initial 
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data for black hole binaries with linear momenta along 
the z-axis, i.e. corresponding to a head-on collision. Even 
though the number of grid points required for reaching 
a given threshold accuracy increases with Z?, we observe 
rapid convergence for all values D — 4, 5, 6, 7. Although 
for large resolutions n the convergence becomes algebraic 
due to logarithmic terms in the r dependency of the so- 
lution u for D = A and D = b, this transition can be 
compensated with a moderate increase in the number of 
grid points n as has also been observed in Ref. [l| for 
Z? = 4. Closer investigation of the profile of the function 
u thus obtained confirms the expected higher falloff rate 
as r — >■ oo for larger D. We further note that u shows 
smooth behaviour near the puncture. Finally, we have 
studied the Hamiltonian constraint as a function along 
the collision axis. As in Z? = 4, the residual constraint 
violations after the elliptic solving are largest near the 
puncture and rapidly falloff away from the puncture. As 
expected from the asymptotic behaviour of u, the con- 
straint violations decay even faster away from the punc- 
ture as D increases. 



The construction of initial data forms a crucial step in 
performing high-energy collisions of black hole binaries 
in higher-dimensional spacetimes which will complement 
existing studies in £> = 4 dimensions [T9'-'2l'] as well as 
studies in Z? = 5 dimensions [18] starting from super- 
posed single black hole initial data. 
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